Dynamic depletion in a Bose condensate via a sudden increase of the scattering length 
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We examine the time-dependent quantum depletion of a 
trapped Bose condensate arising from a rapid increase of the 
scattering length. Our solution indicates that a significant 
buildup of incoherent atoms can occur within a characteris- 
tic time short compared with the harmonic trap period. We 
discuss how the depletion density and the characteristic time 
depend on the physical parameters of the condensate. 



PACS numbers: 03.75.Fi, 05.30Jp 



Bose-Einstein condensates (BEC) of atomic gases with 
tunable self-interaction strengths can be realized by ap- 
plying strong magnetic fields near a Feshbach resonance 
[Q,^). Recent experiments have demonstrated dramatic 
dynamic features by manipulating the values of s-wave 
scattering lengths. For examples, a change of the in- 
teraction sign from repulsive to attractive can trigger a 
condensate collapse HQ , a sweep of the magnetic fields 
across a resonance leads to an enhanced rate of inelastic 
collisions @j. More recently Claussen et al. employed a 
magnetic field pulse to turn on a large scattering length 
in an extreme short time ||. Their observation of con- 
densate loss and its counterintuitive dependence on the 
pulse length suggest rich microscopic dynamics not cap- 
tured by the mean field theory. 

In this paper we address the interesting question how a 
sudden increase of the self-interaction strength affects the 
macroscopic quantum coherence. In particular we will 
examine the dynamic creation of incoherent atoms from 
the condensate (i.e., a time-dependent depletion effect) 
based on a two-body interaction model. Unlike the study 
of stationary systems , the calculation of depletion in 
nonstationary systems is a more difficult task because 
one needs to follow the time evolution of all collective 
excitation modes (assuming these modes are well defined 
in the remote past when the system is stationary) || . In 
this paper we present an approximate solution to tackle 
a time-varying scattering length problem. We find that 
the incoherent atom density can quickly build up if the 
change of scattering length is sufficiently large and rapid. 
Therefore a sudden change of interaction strength could 
lead to observable decoherence effects in the condensate. 
Our solution indicates the time scale of the depletion 
dynamics and the conditions that determine when the 
mean field theory breaks down. 

To begin we consider a Bose-Einstein condensate of N 



atoms confined in a spherical harmonic potential with a 
trap frequency lot- The time-dependence of the s-wave 
scattering length a(t) is controlled by an external mag- 
netic field pulse. We assume that the magnetic field does 
not cross any resonance, and possible molecule forma- 
tion processes will not be considered. Under the dilute- 
ness assumption of the system na(t) 3 <C 1 (n is the atom 
number density), the interaction between atoms is mod- 
elled by a two-body S— potential. The second quantized 
time-dependent Hamiltonian of our model is given by, 



(1) 



Here i& is the atomic field operator, m is the atomic mass 
and ho = ~|^;V 2 + corresponds to the trap 

Hamiltonian for a single particle. We note that the va- 
lidity of replacing the real interatomic potential by the 
S— potential in time dependent situations is a subtle is- 
sue. However, it is plausible that the S— potential is valid 
if the rate of change of the system is small compare with 
a characteristic rate r s _1 that atoms need to adjust their 
relative wave functions The r s is typically a short 
time depending on the scattering length, we will discuss 
this issue later in this paper. 

In this paper we assume that the condensate is station- 
ary initially with a zero scattering length, i.e., a(0) = 0. 
Such an initial condition can be achieved by turning off 
the self-interaction of the condensate adiabatically via a 
magnetic field [||. At later time, a(t) is quickly raised to 
a peak value a > 0. The rise time interval is sufficiently 
short (but longer than r s ) such that the evolution of the 
condensate is negligible during the turn on. This allows 
us to approximate a(t) in the form 



a(t) w a6(t), 



(2) 



where 9 is the step function. 

We decompose the field operator into a dominant co- 
herent part and a small fluctuation part, i.e., ^ = 
\/N&o + We require the condensate wave function $o 
to obey the time dependent mean field equation, 



ih 



dt 



Jx>$o + ff(t)|*of $o, 



(3) 



where g(t) = Airh 2 a{t)N/m. Then according to the 
Heisenberg equation of the fluctuation -0 obeys the 
equation, 



>2,/,t 



(4) 
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Note that in writing Eq. (4) we have assumed that ip 
is small so that second and higher order terms of ip are 
discarded. 

Our task is to solve the time evolution of the opera- 
tor tp and examine the growth of incoherent atom density 
(jp^(x, t)ip(x, t)^. First we need to specify the initial con- 
ditions for $ an d 4>- Since the system is noninteracting 
until t — 0, the initial condensate wave function $o is 
just the ground state of the harmonic trap. The corre- 
sponding quantum fluctuations at zero temperature is de- 
scribed by a vacuum state \vac) associated with ip(x,0), 
i.e., x/j(x,0)\vac) = 0. Using a standard expansion in 
terms of the single particle eigenfunctions of the trap, 
the initial conditions are given by, 



$o(x, t) w f (x)e 



-i\(x)t 



(12) 



$o(j?,0) = /o(£) 

^(#.0) = ^2b n fn(x) 



(5) 
(6) 



Here b n is the annihilation operator associated with the 
trap eigenfunction /„ obtained from hof n — fuv n fn, with 
/o being the lowest state. The initial state \vac) is the 
vacuum of b n . 

The linear structure of Eq. (4) suggests a general so- 
lution in the form 



t) = Y^ bnU n (x, t) + bly* n (x, t), 



(7) 



where the b n and b^ are time-independent operators de- 
fined in Eq. (6). A direct substitution of (7) in (4) gives 
the equation of motion of the mode functions u n (x, t) and 

v n (x, t): 



in 



^={h a + 2g{t)\^{x,t)\ 2 )u n 



-g{t)$l(x, t)v n 



-g(t)$f(x, t)u n . 
To match the initial condition (6), we set 

U n (x,0) = f n (x) 

v n (x,0) = 0. 



(8) 
(9) 



(10) 
(11) 



Thus the operator equation (4) with the initial condition 
(6) can now be replaced by the c— number equations (8- 
11). For later purpose, it is convenient to introduce the 
quantity X(x) = AitnaN /o(:r)/m that characterizes the 
position dependent self-interaction in frequency unit. 

In this paper we present a solution in the short time 
domain < t < T in which T is defined by T = 1/A(0) 
with WtTo -C 1. In this domain the condensate wave 
function does not have enough time to evolve spatially. 
The condensate can only pick up a self-modulated phase 
arising from the self- interaction, i.e., 



Note that the trival phase factor e - luJ Tt/2 ^ ue ^ ^ e trap 
ground state energy is ignored since wt* <C 1 under the 
stated conditions above. We can easily verify the solu- 
tion by a substitution in Eq. (3). The approximation we 
made is the discard of terms involving the spatial gradi- 
ent of e _lA ^^ . These terms would describe the expansion 
of the condensate but they only become significant in a 
longer time scale « 1/ y / WTA(0) > To, which is beyond 
our specified time domain. We have also performed nu- 
merical tests that show Eq. (12) has good agreement 
with the exact numerical solution of Eq. (3) within the 
specified time domain. 

We can solve Eq. (8) and (9) in a similar fashion. We 
write 



u n (x,t) = f n (x)U n (x,t) + 5 nu (x,t) 
v n (x, t) = f n (x)V n (x, t) + 6 nv (x, t) 

where the functions U n and V n arc defined by, 
tk lt = ( h ^ + 2g(t)\^(x,t)\ 2 )u n 



(13) 
(14) 



dt 



+g{t)$l{x,t)V n 
-g(t)% 2 (x,t)U n 



(15) 



(16) 



subjected to the condition U n (x, 0) = 1 and V n (x, 0) = 0. 
We shall show that the 5 nu and 6 nv are small correction 
terms in the short time domain that we are concerned. 

Let us write down the explicit form of U n and V n . Since 
Eq. (15) and (16) do not involve any gradient operator, 
different spatial coordinates are decoupled. At each posi- 
tion the equations behave as a two-level system that can 
be exactly solved: 



U n (x,t) = e- lX ^ 



r> f-\4. .U n + \(x) . 
cos i l n {xjt — i — — sin \l n {x)t 



n n (x) 



sinQ n (x)t 



(17) 
(18) 



which satisfy the normalization \U n \ 2 — \V n \ 2 = 1, and 
the quantity 



fl n (x) = yJuJn (uJ n + 2A(x)) 



(19) 



can be interpreted as a form of frequency analogous to 
the Rabi oscillations in two-level systems. 

To show S nu and 8 nv are small, we note that they are 
both zero initially. Their equations of motion can be 
obtained by a direct substitution of Eq. (13) and (14) 
into (8) and (9). This leads to a set of linear equations 
similar to Eq. (8) and (9) but with inhomogencous terms 
as 



2 



at m \ 2 

at m \ 2 



(20) 
(21) 



These inhomogeneous terms involve the spatial gradient 
of /„, U n and V n . Since U n and V n are smooth func- 
tions in the length scale of the trap ground state, the 
dominant contribution comes from V/ n for high excited 
modes. We can estimate the size of S nv using the first 
order approximation which is a time integration of the 
inhomogeneous term W n ■ V/„. Up to the time t = Tq, 
we find that S nv is only a small fraction of f n V n . The 



fraction is typically of the order (utTq) 1 ^ 2 -C 1 fl 
similar conclusion applies to S nu . 

The functions f n U n and f n V n are approximate solu- 
tions of u n and v n . The approximation comes from the 
fact that V 2 (f n U n ) pb £/ n V 2 /„ and V 2 (/„K) « KV 2 /„ 
are the leading contributions in the Laplacians for slowly 
varying (spatially) functions U n and V n . The correction 
terms can remain small under the conditions stated in 
the previous paragraph. We have tested the approximate 
scheme beyond the first order consideration by solving 
Eq. (8) and (9) numerically in a one-dimensional model, 
the exact results are in good agreement with those ob- 
tained by the approximate method. 

Now we come to the main result of our paper. The 
time-development of incoherent atom density after the 
change of the scattering length is given by, 



A 



ft (x,t)i>{x,t)j =Y / K(s,t)\ 2 

n 

^^2\fn(S)\ 2 \V n (x,t)\ 
n 

With the help of Eq. (18), we have, 

ft (f,t)^(a,t))«X)LM 



(22) 



u?y 2 



A 2 (f) 



x sin 



LO n (uj n + 2 A (£)) 

\^Juj n (uj n + 2A (x))t 



(23) 



For a spherical harmonic trap, each single particle state 
n is described by three quantum numbers (nj,, n^, TI3) in 
Cartesian coordinates, and the corresponding frequency 
w n = (ni + n 2 + n 3 + 2>/2)lot- Therefore Eq. (23) is a 
triple summation. 

Let us calculate the incoherent atom density at the 
trap center x — 0. It is a good approximation (particu- 
larly for high modes) that 



|/2ji,2j 2 ,2j 3 (0)1" 



2ji 



0.09 x I- 3 



2h + 5 



2,3 + - 2 



1-1/2 



(24) 



where Iq — y/h/rrujr and we need only to consider even 
modes at the origin. Eq. (24) is obtained from the semi- 
classical construction of the trap eigenfunctions. One 



can check that Eq. (24) is fairly accurate even for the 
ground state. Next we change the variables ji by in- 
troducing: k 2 = (4ji + l)w T T /2 for i = 1,2,3. We 
find that the function inside the summation (23) changes 
smoothly with fej. Therefore we may replace the dis- 
crete sum (23) by an integral that can further be simpli- 
fied using the spherical symmetry. In order to compare 
the fluctuations with the mean-field density, we define 

A(t) = ^^(0, t)V>(0, t)J /no which is the ratio of the cen- 
tral incoherent atom density to the condensate density 
n Q = iV/ 2 (0). We find that 



A(t) « 6.3(n a 3 ) 1/2 / dk 



sin 2 (kty/2 + k 2 /T 



k 2 



(25) 



where k 2 = kf + fc| + k 2 . The integral in Eq. (25) is 
an increasing function of time, but it must be less than 
7i"/2v2 w 1.1. At the time t — Tp the integral reaches 0.53 
that gives A pa 3.4(n a 3 ) 1/2 JlJ. This result is about 
twice the value that one would obtain from a (station- 
ary) homogeneous condensate at zero temperature [^2[ . 
Note that A(t) depends explicitly on the gas parameter 
noa 3 in Eq. (25). If A(i) is comparable to one, then both 
the mean field theory and the linearized quantum theory 
break down because quantum fluctuations are no longer 
small. For a general position x, we can work out a sim- 
ilar integral. The spatial profile can be well fitted by a 

Gaussian: lft{x,t)<tp(x,t)\ pa (ft(0,t)fto,t)\ e" 3r2 / 2 'o. 

Finally we would like to discuss the validity of the 
model (1). A necessary condition that the usual 
pseudopotential approach can be generalized to time- 
dependent situations is that the time- varying parameters 
of the system do not excite high relative wave vectors be- 
tween atoms. This condition ensures pa/h <C 1 (where 
p is the relative momentum) required for the s-wave ap- 
proximation. There is no simple way to check the con- 
dition for general interatomic potentials, but it is useful 
to consider the hard sphere model as an illustrative ex- 
ample. In the relative coordinate of two hard spheres, 
the expanding spheres are described by a moving bound- 
ary condition. If the radial velocity is d, a slow incom- 
ing wave will pick up a wavevector of magnitude roughly 
m r a/h by reflection, where m r — m/2 is the reduced 
mass. Therefore m r aa/h <C 1 is the required condition 
that specifies how slow the speed should be. If we con- 
sider a(t) = at/trise is linear in time during a rise period 
t r i se , then the condition is satisfied if m r a 2 /h <C t r i se . 
Quite remarkably, the characteristic time t s = m r a 2 /h 
estimated here agrees with the one suggested by Leggett 
under a general consideration ||. 

To give a physical example, we consider a 85 Rb conden- 
sate with 10 5 atoms in a spherical trap with lot — 27r x 10 
Hz. If the scattering length a = 30 nm, then To = 80/xs 
and t s = 1/is is a short time. In this example, the gas 
parameter n a 3 = 0.001 and A = 0.11 at t = T . This 
means a 11% depletion of condensate density at the trap 



3 



center in 80/is. We may compare this example with the 
case with a smaller scattering length a — 3 nm (same par- 
ticle number and trap frequency) that only gives 0.36% 
depletion density in 800/its. 

In conclusion we made use of the model (1) to ana- 
lyze the dynamic depletion in a Bose condensate after 
a sudden increase of the self-interaction strength. The 
depletion is a quantum decoherence process originated 
from the vacuum fluctuations associated with ip. Our 
solution (23) and (25) indicate how the depletion den- 
sity increases as a function of time. We discover that 
To = m I AnhaN f§ (0) is an important time scale for the 
buildup of incoherent atoms for a trapped condensate 
with ujtTo C 1. It is interesting to rewrite To in the 
form: T = Voa 3 )~ 1/3 (™0© _2/3 / 47r - This indicates 
the relationship between the characteristic time To and 
different parameters, including the gas parameter. Fi- 
nally, we point out that the recent experiment at JILA 
H reaching noa 3 « 0.3 is beyond the scope of the present 
theory. Such a strong coupling would give t s w 39/zs that 
is longer than the rise time of the magnetic pulse, and so 
the s-wave approximation model is questionable. The mi- 
croscopic processes in such an extreme situation remain 
an interesting open topic for investigations. 



A(2T ) » 3.5(n a 3 ) 1/2 which is only about a 3% increase 
compared with the value at t = To. 
[12] K. Huang, Statistical Mechanics (Wiley, New York 1987). 
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